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Abstract 

The electromagnetic field inside a cubic cavity filled up with a lin- 
ear magnetodielectric medium and in the presence of external charges 
is quantized by modelling the magnetodielectric medium with two 
independent quantum fields. Electric and magnetic polarization den- 
sities of the medium are defined in terms of the ladder operators of 
the medium and eigenmodes of the cavity. Maxwell and constitutive 
equations of the medium together with the equation of motion of the 
charged particles have been obtained from the Heisenberg equations 
using a minimal coupling scheme. Spontaneous emission of a two level 
atom embedded in a magnetodielectric medium is calculated in terms 
of electric and magnetic susceptibilities of the medium and the Green 
function of the cubic cavity as an application of the model. 

Keywords: Field quantization, Magnetodielectric, Spon- 
taneous emission, external charges 

PACS number: 12.20.Ds 

1 Introduction 

The usual approach to quantum optics in non dispersive dielectric media is to 
introduce the medium by its linear or nonlinear susceptibilities. In one proce- 
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dure, the macroscopic fields are used to build an effective Lagrangian density 
whose Euler-Lagrange equations are identical to the macroscopic Maxwell 
equations in the dielectric, where the constitutive equations couple the po- 
larization field to the electric field However, attempts to add disper- 
sion to this effective scheme have run into difficulties [7[ E]. The reason for 
this is that inclusion of dispersion leads to a temporally nonlocal relation- 
ship between the electric field and the displacement field and the effective 
Lagrangian in this case is also nonlocal in time and cannot be used directly 
in a quantization scheme. Another problem is the inclusion of losses into 
the system. It is well known that the dissipative nature of a medium is an 
immediate consequence of its dispersive character and vice versa according 
to the Kramers- Kronig relations [S]. This suggests that in order to quantize 
the electromagnetic field in a dielectric in a way that is consistent with the 
Kramers-Kronig relations, one has to introduce the medium into the for- 
malism explicitly. This should be done in such a away that the interaction 
between light and matter will generate both dispersion and damping of the 
light field. In an attempt to overcome these problems, by taking the po- 
larization of the medium as a dissipative quantum system and based on the 
Hopfield model of a dielectric [TOlfTT] . a canonical quantization of the electro- 
magnetic field inside a dispersive and absorptive dielectric can be presented, 
where the polarization of the dielectric is modelled by a collection of interact- 
ing matter fields [T^ITB] . The absorptive character of the medium is modelled 
through the interaction of the matter fields with a reservoir consisting of a 
continuum of the Klein-Gordon fields. In this model, eigen-operators for cou- 
pled systems are calculated and the electromagnetic field has been expressed 
in terms of them. The dielectric function is derived in terms of a coupling 
function which couples the polarization of the medium to the reservoir. 

One approach to study a quantum dissipative system tries to relate dis- 
sipation to the interaction between the system and a heat bath containing 
a collection of harmonic oscillators [ll]-[26]. In this method, the whole sys- 
tem is composed of two parts, the main system and a bath which interacts 
with the main system and causes the dissipation of energy on it. In order to 
quantize the electromagnetic field in the presence of an absorptive magne- 
todielectric, it is reasonable to taking the electromagnetic field as the main 
quantum dissipative system and the medium as the heat bath. In this point 
of view the polarizability of the medium is a result of the properties of the 
heat bath and accordingly, the polarizability should be defined in terms of 
the dynamical variables of the medium. In this method, contrary to the 
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damped polarization model, polarizability and absorptivity of the medium 
are not independent of each other [121 [IS]- Furthermore, if the medium is 
magnetizable, as well as polarizable, the medium can be modelled with two 
independent collection of harmonic oscillators, so that one of these collec- 
tions describes the electric properties and the other describes the magnetic 
properties of the medium. This scheme propose a consistent quantization of 
the electromagnetic field in the presence of an absorptive magnetodielectric 

In the present paper, we generalize the quantization scheme presented in 
[27] to the case where there are some external charges in the medium. We 
take a cubic cavity filled up with a linear magnetodielectric, containing some 
external charges, as the medium. As a simple application of the method, 
we study the spontaneous decay of an initially excited two-level atom in a 
dispersing and absorbing magnetodielectric. 



2 Quantum dynamics 

Quantum electrodynamics in a linear magnetodielectric medium, can be ac- 
complished by modelling the medium with two independent quantum fields 
that interact with the electromagnetic field. One of these quantum fields, 
namely E-quantum field, describes the polarizability character of the medium 
and interacts with displacement field through a minimal coupling term. The 
other quantum field, M-quantum field, describes magnetizability character 
of the medium and interacts with magnetic field through a dipole interaction 
term. The Heisenberg equations for electromagnetic field and the E- and M- 
quantum fields, give not only the Maxwell equations but also the constitutive 
equations of the medium and the equations of motion of external charges. 

For this purpose let us consider an ideal rectangular cavity with sides 
Li, L2 and L3, along the coordinate axis x, y and z, respectively, filled up 
with a magnetodielectric. The vector potential of the electromagnetic field 
in coulomb gauge can be expanded as [28j 
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where 
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h{n,r) = sm— — sm — — cos— — , (2) 
Li L2 -L3 

the mode vector n is a triplet of natmal numbers (ni, n2, 77-3), means 
n2 n3=i' ^ ~ -^1-^2-^3 IS the volume of the cavity, Eq is the permittivity 

of the vacuum, ujfi = c^J + + ^^||- is the frequency corresponding to 
the mode n and e(n. A), (A = 1, 2) are polarization vectors which satisfy 

e{n, X) ■ e{n, X') = ^^a', 

e(n. A) • = 0, (3) 

where A;^ = + ^j^j + ^^A; is the wave vector. Operators a^x{t) and 

^Hxi^) annihilation and creation operators of the electromagnetic field 
and satisfy the following equal time commutation rules 

[anx (t) , aiy (t)] = SH,mSx\' ■ (4) 

— * — * 

Let F{r,t) be an arbitrary vector field, the transverse component F^{f,t) 
and the longitudinal component Fll(f, t) of F{r^t) are defined as 

F^{r, t) = F(f, t)+ [ dh'V ■ F{P, t)VG{f, r'), (5) 
F\f,t) = - j c/VV • F{r',t)VG{f,r'), (6) 

where G{f,r') is the Green function of the cubic cavity 

-ix 8^1 . niTvx , n2'ny . nsnz . riinx' . n2'ny' . n^nz' 
G[r, f ) — T7 — — sm — — sm ^ — sm — — sm — ; — sm — ; — sm ■ 



^ n l^nP -^3 -^1 -^2 -^3 

(7) 

In the presence of external charges the displacement field is not purely trans- 
verse. The transverse component of the displacement field can be expanded 
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in terms of the cavity modes as 
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n A=l 



The commutation relations between the components of the vector potential A 
and the transverse components of the displacement field, D-*-, can be obtained 
from (HI) as follows 



[Ai{r,t),-Df{r',t)] = thS^^{f-r'), (9) 
t 

eigenvector fields m^a 



where 5[^fr — r') is the transverse delta function defined in terms of the 



n A=l ft \'^n\ 

(10) 

The Hamiltonian of the electromagnetic field inside a magnetodielectric in 
the normal ordering form is 

2 

= im^n4AW«nA(t), (H) 
n A=l 

where /xq is the magnetic permittivity of the vacuum and : : denotes the 
normal ordering operator. By modelling the magnetodielectric medium with 
15- and M-quantum fields we can write the Hamiltonian of the medium as 
the sum of the Hamiltonians of the E- and M-quantum fields 



Hd = Hp + H 



m 

oo 



d^kfkj^,dUk,t)driAk,t) 
n v=i 

3 c+oo _, _, 

Hmit) = Y.T. d^knw^bi,ik,t)brAk,t), (12) 

1 J —oo 



n u=l ' 



where a;^ is dispersion relation of the medium. Quantum Dynamics of a dis- 
sipative harmonic oscillator interacting with an absorptive environment can 
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be investigated by modeling the environment with a continuum of harmonic 
oscillators [I1|-[IH]- In the case of quantization of electromagnetic field in the 
presence of a magnetodielectric, electromagnetic field is the main dissipative 
system and the medium play the role of the absorptive environment. The 
Hamiltonian ffTTl) shows that electromagnetic field contains a numerable set 
of harmonic oscillators labeled hj n, A. Therefore to each harmonic oscillator 
of the electromagnetic field labeled by n, A, a continuum of oscillators should 
be corresponded. In the present scheme, to each harmonic oscillator of the 
electromagnetic field labeled by n, A, we have corresponded two continuous 
set of harmonic oscillators defined by the ladder operators dfiu{k, t), dl^^^k, t) 
and bfi^{k,t), bl^{k,t) which are to describe the electric and magnetic prop- 
erties of the medium respectively and satisfy the equal time commutation 
relations 

[dr{u{k,t),dl^,{k',t)] = Sri,mSuu'S{k - k') , 

[bHu{k,t),bl^,{k',t)] = 6H,rnSuu'S{k-k'). (13) 

If we want to extend the Huttner and Barnett model [T2| to a magnetodi- 
electric medium, we have to take a vector field X as electric polarization 
density of the medium and a heat bath B interacting with X in order to take 
into account the absorption due to the dispersion induced from the electrical 
properties of the medium. Similarly, we must take a vector field Y as mag- 
netic polarization density and a heat bath B{ independent of B) interacting 
with Y in order to take into account the absorption due to the dispersion 
induced from the magnetic properties of the medium. Following the Huttner 
and Barnett method, we find that a Fano diagonalization [TT] process of the 
sum of the Hamiltonians related to the electric polarization X, the heat bath 
B and their interaction, lead to what we have named H^, in Eq. f[T^ . Also a 
Fano diagonalization process of the sum of the Hamiltonians related to the 
magnetic polarization Y, the heat bath B and their interaction, lead to what 
we have named Hm in Eq. f|T2|) . If we use the Huttner and Barnett approach 
inside a rectangular cavity, the labels n and u will appear in the expansion 
of the electric polarization field X in terms of the eigenmodes of the cavity. 
On the other hand, the heat bath B interacting with the polarization field X 
consists of a continuous set of harmonic oscillators labelled by a continuous 
parameter u. It should be noted that the three labels n, v and u; remain in 
Bf, after the diagonalization process. 

For a linear medium, the electric polarization density operator can be 
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written as a linear combination of ladder operators d^y{k^ t) and d\^{k, t) 

P{r, = y ^ E E / d'k[f{uj;, r)driu{k, t) + h.c.]vrUr) (14) 



ri !^=1 

with 



VHu{.r) = Uriu{r), z/ = l,2 

k 



Vnsir) = k^^fi{n, f)x + A;^y/2(n, f)y + k^Jsin, r}z, kft 



(15) 



the function /(w^, r) is the coupling function of electromanetic field and the 
E-quantum field. Also for a linearly magnetizable medium we can express the 
magnetic polarization density operator of the medium as a linear combination 
of the ladder operators of the M-quantum field [27]. 

M{f,t) = w|rEE / d'k[g{u-„r)briu{k,t) + h.c.]sriu{f^. (16) 

where the function g{ujj^, r) couples the electromagnetic field to the M-quantum 
field and 

Snv{r) = ^ ^ ^"^ , 1^=1,2, 

\kn\ 

SH3ir) = knxQiin, r)x + hy92iri, r)y + knzQsin, r)z, 
gi {n, r) = sin ■ 

92 {n, r) = cos 
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g3{n,-r) = cos— — cos— — sin— — . (17) 
Li L2 L3 

Now let the total Hamiltonian, i.e., electromagnetic field , the 15- and M- 
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quantum field and the external charges be like this 

H{t) = (frl ^ ^ ' + ^ ' ^ ' ' - V X A((f, t) ■ M{f, t) 

Jv I 2^0 2/io 

+ H„, + }^ + 7^1^ qaqpG[r^, rp) 

N 



— E^- / rfVG(f„,rOV'-P(r',t) 



where g^, m^, and Tq, are charge, mass, linear momentum and position of 
the ath particle respectively. The function G is the Green function given by 
(I?!). Using ([9]), the Heisenberg equation for A and can be obtained as 

g^4l^.i(r.«)|^- ^^<--'> (19) 
at n fxo 

(20) 

where 

^ ^ (if 

JHr,t) = T.T.Qo.-^^Hirc.{t) - r) (21) 

is the transverse component of external current density. If we define trans- 
verse electrical field i?"*", induction B and magnetic field H as 

dA B 

E^ = - — , B = VxA, H= M, (22) 

dt fio 

then (|T9|) and (|20|) can be rewritten as 

= £o^^ + P^, (23) 

__ + = V X if, (24) 

as expected. In the presence of external charges, the longitudinal components 
of the electrical and displacement field can be written as 

1 ^ ^ Pll 
E\\{f,t) = E gaV.G(f,r;(t)) , (25) 
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N 



D\r,t) = e^eW +p\\ = -Y, qa^rG{f,Ut)), (26) 

a=l 

where P" is the longitudinal component of electric polarization density. If 
we apply the Heisenberg equation to the ladder operators of the medium, we 
obtain the macroscopic constitutive equations of the medium which relate the 
electric and magnetic polarization densities to electric and magnetic field, 
respectively. Using commutation relations fll3l) . we find from Heisenberg 
equations 

dfiu{k,t) = - [H,dH^{k,t)] = 

^1^1 ,3 J f*. 



-lu-^dnAk.t) + J^d'r'r{uj-^,r')E{r',t) ■ t7^,(r'), (27) 



hHv{k,t) = -[H,bjtu{k,t)] = 
- lu^hnuik, t) + 1^ J^d'r'g*{u;^, r')B{r', t) ■ sn^i?), (28) 

with the following formal solutions 

rf^,(^,t) =rf^,(fc,0)e-^"£* + 

f dt'e-'->^^'-''^ I d'r'nu^,, r')E{r', t') ■ vU?), (29) 
h\ V Jo Jv 



^/f ll jy d'^'9*i^k^ ^')B{r', t') ■ sn^i?). (30) 

By substituting (129|) in (IT^ and using the completeness relations 

E E v'^Ar-)vi^{r') = E E ^;^.(^14.(^0 = jS^^Sif- ?), (31) 

fa u=l in v=\ 

we can find the constitutive equation 

P(f, t) = P^if, t) + Bo /'*' dt'xeir, \t\ - t')E{f, ±t'), (32) 

Jo 
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for the medium which relates the electric polarization density to the total 
electric field, where the upper (lower) sign, corresponds to t > (t < 0) 
respectively, and E = E-^ + -E" is the total electrical field with E-^ and 
defined in f l22]) and fl2^ . The memory function 

Xeir,t) = { (33) 

t < 

is the electric susceptibility and in frequency domain satisfies the following 
relations 

Im[x {r,u) = —\fiu,f)\^y2- L'' , (34) 

Re[x^{f,u) = — d\k\\kmu;^„r)\'^^^ (35) 



where 



X^{f,u)= I dtxe{r,t)e'''' (36) 







is the Fourier transform of the electric susceptibility and \ ki\s are the roots of 
the algebraic equation u = = u{\k\). A feature of the present quantization 
method is it's fiexibility in choosing a dispersion relation and a coupling func- 
tion such that they satisfy the relation flMl) . In other words, the dispersion 
relation and the coupling function /(cug, r) are two free parameters of this 
theory up to the constraint relation ( IMI) . Various choices of a;^ and /(o;^, r) 
satisfying do not change the commutation relations between electromag- 
netic field operators and lead to equivalent expressions for the field operators 
and physical observables. As can be seen from the relation in order to 
have a finite value for susceptibility for any frequency a;, we should assume 
that the denominator of the fraction in is non zero, so we assume that 
the dispersion relation is a monotonic function of . For simplicity and 
simple calculations we take a linear dispersion relation uj = c\k\ which leads 
to the following relation for the electric susceptibility: 



Xeir,t) 



f J^/o°°da;a;2|/(^,OPsma;t t>0 

(37) 



t < 
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The operator P^ir^t) in fl32p . is the noise electric polarization density 

Piv(r,t) = Mt.T. [ d'k[f{u^,r)drt,{k,0)e-'^^^' + h.c.]vrtAf^. (38) 



Similarly by substituting fl30l) in f|T6l) . we obtain the following expression for 

M(f, t) 



M{f, t) = M^(f, t) + — f dt'xmir, \t\ - t')B{f, ±t'), (39) 

/iO JO 

where Xm is the magnetic susceptibility of the medium 



Xml'",^) = duio \g{uj,r)\ smut, t > 0, 

xUr,t) = 0, t<0. (40) 

and operator Mjs[{r,t) is the noise magnetic polarization density 

Mjv(r,t) = Mt.T.J d'k[g{uj^,,r)bn.{k,0)e~'->^' + H.C]srt,{f^. (41) 

' n 1/=! 

If we are given a definite Xe{t) and Xm(^), then we can inverse the relations 
( |37|) and P3|) and obtain the corresponding coupling functions / and g as 

|/(a;,f)|2 = dtXe{r,t) sin ut, uj > 0, 

An^oj^ Jo 

|/(a;,r-)|2 = 0, u = 0. (42) 



\g{iO,f)\'^ = 4^2^^^2 /q c?tXm(^,^) sinwf, a; > 0, 



Therefore the constitutive equations (]23l) . ( 1321) and (1411) together with the 
Maxwell equations, can be obtained directly from the Heisenberg equations 
applied to the electromagnetic field and the quantum fields E and M. 
It is clear from (|38|) and (HTj) that the explicit forms of the noise polarization 
densities are known. Also because the coupling functions /, g are common 
factors in the noise densities Pn,Mn, and susceptibilities Xe,Xm, it is clear 
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that the strength of the noise fields are dependent on the strength of Xe^ Xmi 
which describe the dissipative character of a magnetodielectric medium. Spe- 
cially when the medium tends to a non absorbing one, the noise polarizations 
tend to zero and this quantization scheme is reduced to the usual quantiza- 
tion in this medium [27]. 

Finally using the total Hamiltoian flTSl) and the commutation relations 

[ra.Pp] = iMapI (44) 

it is easy to show that the Heisenberg equations of motion for charged par- 
ticles are 

rriafa = qaE{fa,t) + 

^[f, X B{f^, t) + B{r^, t) X f,] a = 1, 2, iV (45) 

3 Spontaneous emission of an excited two- 
level atom in the presence of a magnetodi- 
electric medium 



In this section as a simple application we use the quantization scheme in 
the previous section to calculate the spontaneous decay rate of an initially 
excited two-level atom embedded in a magnetodielectric. For this purpose, 
let us consider a one electron atom with position R which interacts with the 
electromagnetic field in the presence of a magnetodielectric medium. If we 
restrict our attention to the electric-dipole approximation [29l 130], then the 
Hamiltonian flTSl) can be approximated as 



H = Ho + H\ 

Hq = Hp + He + Hrn + Hat, 



H' 



V 



V xA- M 



ef-ED{R,t), 



(46) 
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where 



-2 ^ 

Hat = \ qaG{r,fa), 



t) = () - . _ f(M, (47) 

£o £o ^0 

and e, m, r and p are charge, mass, position and momentum of the atom elec- 
tron. In Hamihonian (l46l) we have ignored from j^A^{R, t) and 2^ d^rP'^{f, t) 
, because these terms do not affect the decay rate or level shifts in the dipole 
approximation and therefore will not be of concern here. For reality we as- 
sume that the atom is localized in a free space region Vq without polarization 
density. Then for a two-level atom with upper state |2), lower state |1) and 
transition frequency ujq, the Hamiltonian can be written as ^28j 

H = hujoa^a + Hp + + + Hat 
■ P 



3^' 



+ / d^r 



^0 



-V X A-M 



^0 



(48) 



where Q = V — Vq, is the region containing the magnetodielctric, a = |1)(2| 
and cr^ = |2)(1| are Pauli operators of the two-level atom and fi2 = {^\r\2) 
is its transition dipole momentum. It is remarkable that the expansions (HM 
) and flTBl) for electric and magnetic polarization densities are independent 
of the shape of the magnetodielectric and depend only on the shape of the 
cavity contained the magnetodielectric. The shape of the medium affects 
the coupling functions f{uj,r) and g{uj,r) and the time dependence of the 
ladder operators of the medium, such that for a given medium with definite 
susceptibilities Xe and Xm, the dynamics of the total system leads to the 
correct constitutive equations fl32l) and fl39|) where the coupling functions are 
given by fj42|) and fH3l) . This is just as the expansion ([1]) for the vector 
potential of the vacuum field where the presence or absence of the external 
charges do not change the form of the expansion and the presence of the 
external charges merely affect the time dependence of the ladder operators 

To study the spontaneous decay of an initially excited atom we may look 
for the wave function of the total system in the framework of the Weisskopf- 
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Wigner theory [ 

m)) = c{t)\2)\0)F\0),\0U + J2j2Fnxm)\n,X)F\0)e\0)r 

n A=l 

+ EE J d''kD^,ik,t)\l)\0)F\m,k,u),\0)^ 

m u=l 

+ EE J d^kM,r.Ak,t)\imF\0)e\m,k,u)r^, 



(49) 



where |0)ir, |0)e and |0)m are vacuum state of electromagnetic field and E- 
and M-quantum fields, respectively. The coefficients c(t), F^a(^), D^^{k,t) 
and M^j,{k,t) are to be specified by Schrodinger equation 



(50) 



for initial conditions c(0) = 1,F^a(0) = Drfw{k,Q) = M^,^(A;, 0) = 0. Substi- 
tuting I'lpit)) from ( H9l) in ( l50l) and using the expansions ([H]), (JHI), ( |T6l) and 
applying the rotating wave approximation(RWA) [29l [30], we find the fol- 
lowing coupled differential equations for the coefficients of the wave function 



n A=l 

ihb^,{k,t) = nw^,D^,{k,t) + Y.il[QtxH)TFnx{t) + 

n A=l 



(51) 



(52) 



m v=l 

EE fd'kQ^a^k)DrnAk,t) + 



m 1^=1 



c{t), 



ihc{t) = hjjQcit) 



%e 



Fnx{t), 



(53) 



(54) 
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where 



d^r'f{uj^, r')uH\{r') ■ v,^u{r') 



1 32h^oUJfi 



Jv-Vo 



d r'g{ujj:,r')sfix{r') ■ SrAu{r'). 



(55) 



One can solve these coupled differential equations by Laplace transformation. 
Let f{p) denotes the Laplace transformation of f{t), then taking the the 
Laplace transformation of equations fl5Tl) - fl53l) and then their combination 
and using the completeness relations fl3T|) we find 



„-/ A'=l 



■ri2 ■ uh\{R) 



(56) 



where we have applied the initial conditions F^a(O) = D^fiuik, 0) = Mj^^^k, 0) 
and 



Vy ujf, 



X 



d r 



Ze{ip,r)un\{r) ■ u-,y{r) + Zm{ip,r)sH\{r) ■ s^,y{r) 



(57) 



Ze{ip,r) = ^ [ du ^'''^^'^\ Zm{ip,r) = I du^""'^^'^^ 



271 Jo ip — U 271 Jo ip — UJ 

(58) 

Here Xeii.^, and Xmi(<^, "T) are the imaginary parts of electric and magnetic 
susceptibilities in frequency domain, respectively. Equation (15^ is a compli- 
cated algebraic equation for F^a(p) and may be solved by iteration method. 
In the first step of an iteration method may be approximated by 

FM - (P) + E E ,^_,_,|i + M,2.,,p„. (59) 

15 



where 



le 



-c(p), 



(60) 



is the solution of (156|) in the absence of medium, that is when W^';^' = 0. 
Now combination of equation (159!) and the Laplace transformations of the 
equation ( 1541) yields 



ih{c{p) - c(0)) = nwoc{p) + [r[^ ■ {G^°\R, R, ip) + G{R, R, ip) 



■ ri2 



where 



c(p), 



(61) 
(62) 



G{R, R, ip) 



4e2 



E E 



^,K,ir„^,(i?)w^;^'(zp)«„,,,(i?) 



A,fei [^p - (i + w^>:{ip) 



ip - 



(63) 



The diadic G'-*'-* gives the contribution of the vacuum field for the spontaneous 
emission and level shift of the atom, that is, when there is no magnetodielec- 
tric, while the diadic G gives the spontaneous and frequency shift of the atom 
due to the presence of the medium. From definitions flFTI) and flSHl) it is obvi- 
ous that Ze(zp, r), Ze{ip,r) and W^'^'^' (tp) are analytic functions in complex 
semi plane Re{p) > and since the denominator in diadic G has no zero in 
semi plane Re{p) > so G{R, R, ip) is analytic for Re{p) > 0. Taking the in- 
verse Laplace transformation of (!60|) we find the following integro-differential 
equation 



c(t) = -lujQcit) - (Fo + 2Ao)c(t) + / Kit- t')c(t')dt', 

Jo 



where 



K{t - t') 



1 

2^ 



fi2- (G{R,R,u + ir)) -fia 



(64) 



(65) 



and Fq and Aq are contribution of the vacuum field for the decay rate and 
the lamb shift of the two-level atom [HT]. Here we restrict our attention to 
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the weak coupling regime, where the Markov approximation j5T| [52] apphes. 
That is to say, we may replace c{t') in the integrand in fl61l) by 

c(t') =c(t)e^""(*-*'), (66) 

and approximate the time integral in (1611) in very large times by 



iP h tt6{ujo — u] 



(67) 



where P denotes the principal cauchy value. After some simple algebra and 
using the Kramers-Kronig relations for the diadic G we deduce 



c{t) = -iuJoc{t) - (Fo + r + lAo + zA)c(t), (68) 



where 



Vi2 ■ ImG{R, R, uo + lO^) ■ fu 



h 

1 



ri2- ReG{R,R,uJo + iQ^)-r 12 , (69) 



are decay constant and the level shift due to the presence of the magnetodi- 
electric, respectively. 

There are many cases that the medium is homogeneous inside the region 
Vt and W^^{ip) in the denominator of G in fl63|) is independent of polarization 
label A and vector label n. For example for the box 0<a;<Li,0<|/<L2, 
< z < -y- filled up uniformly with a homogeneous magnetodielectric, from 
the definition (157|) . we have 

Wi^iip) ^ W{ip) = \ {Z^ip) + Z„(zp)) . (70) 

Generally speaking for a piece wisely homogeneous medium with homo- 
geneous pieces in volumes fii, ^2-,- ■ ■ , ^n, applying the definitions ([2]) and ( 
[T7|) in fl57|) . we note that W^^{zp) in denomonator of G can be approximated 

by 

N r) 

WI>:M ^ W{^p) ^ E y {Zj^\^p) + Z^^^W)) > (71) 

where and Zf^ denote the quantities and Z^ for the ith piece respec- 
tively. Now we can substitute W^'^' {ip) from ( 1571) in the numerator of G 
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given by (1^ and do the summations over A, A' by using the completeness 
relations 

A=l 

2 ^ 3 ^ 

E 4a(^)4a(^) = E ^jt^iknjiin, R)gj{n, r). (72) 
A=i ^J.=l 

Now in the box frame we have 

(73) 

where 

<,{R. r, ^P) = ^T. _ l^"^ ^(,^)] (^«7 - k-,kt) Un, i?) A(n, f), 



C|^(^, f, «p) = E ^" {^Pi - ^^H^t) //3(^> R)fy{n,r), 



4 3 



V^^iR, f, ^p) = E E -^^7M/3^^^//3(^, R)g,in, r-) (74) 

It is remarkable to note that since we have assumed that the center of the 
atom is localized in a free region, then in the argument of the tensors rj^, 
C^, rj"^ and in dZlD, we have R ^ f and therefore the recent series are 
all convergent. Inspection of the definitions (l58ll shows that the real and 
imaginary parts of Ze{uj + lO'^) and ZmioJ + tO'^) for real frequency u are 



= aeiuj) 



1 /""^ 

Re[Ze{uJ + tO+)] = —P / rfa; 

ZTT JO 

/m[Ze(u; + «0+)] = -^Xei(^)^(^) = 7e(t^), 

Zn Jo uj — 00 

Im[Zm{iO + lO'^)] = -^Xmi(^)^9(^) = 7m(w), (75) 

18 



where 'd{u) is the step function and P denotes the Causchy principal value. 
Finally using dnSD , (IZ3D, dZi), dZSD and 



1 



1 



in6{iL! — uj'), 



(76) 



uj + iO+ - uj' 



P 



we find the spontaneous emission and the shift frequency of a two level atom 
due to the presence of a magnetodielectric. From (!69|) and (!73|) it is clear 
that the spontaneous emission and the lamb shift depend on the (i) position 
of the atom, (ii) orientation of the dipole of the atom,(iii) size and sort of 
the medium, (iv) relative position of the medium respect to the center of 
the atom, (v) the volume and the shape of the cavity. For a a non rect- 
angular cavity we must use a different set of eigenvectors and therefore the 
spontaneous emission depends also on the shape of the cavity. 

4 Conclusion 

The scheme introduced by the present authors to quantize the electromag- 
netic field in a magnetodielectric is generalized to the case where some ex- 
ternal charges are present in the medium. The spontaneous emission of a 
two-level atom embedded in a magnetodielectric is calculated as an applica- 
tion of the model. 
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